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Abstract 

In the light of KamLAND data just released, we reanalyze and update the constraints on neutrino 
' masses and mixing parameters, the most general ones that can be derived in three-flavor mixing 

■ scheme of neutrinos with use of the bounds imposed by neutrinoless double beta decay search and 

Psj ■ reactor experiments. We point out that with KamLAND data and assuming Majorana neutrinos 

one can derive, for the first time, an upper bound on neutrino contribution to the cosmological 
parameter, VL^ < 0.070/i~^, by using the current upper bound on mass parameters obtained by 
Heidelberg-Moscow group, {m)pp < 0.35 eV (90 % CL). If the bound is tighten to {m)pp < 0.1 eV 
by future experiments, it would lead to a far severer bound fij, <0.01/i^^. 
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I. INTRODUCTION 



The solar neutrino problem is finally solved by the KamLAND experiment which 
observed a clear deficit of about 40 % in z/g fiux from reactors located in all over Japan. 
It pin-points the Large-Mixing-Angle (LMA) region as the unique solution of the solar 
neutrino problem based on the MSW matter enhanced neutrino fiavor transformation 
excluding the remaining LOW and the vacuum (VAC) solutions. While the LMA solution 
was preferred over others on statistical ground in combined analyses of solar neutrino data 

it is by no means true that it was selected out as the unique solution. It is because the 
values of of the LOW and the VAC solutions, ~ unity per degrees of freedom^, were 
acceptable and hence it was not possible with only solar neutrino data to uniquely idetify 
the LMA solution. 

It implies, first of all, that the three-fiavor mixing scheme of neutrinos had further support 
by the experiments, and the structure of lepton fiavor mixing is established together with 
the existing evidences which come from the atmospheric and the long-baseline accelerator 
experiments [6] . Given our above understanding of the KamLAND result it also means that 
the experiment gave the first direct proof of (almost) pure vacuum neutrino oscillation, the 
fascinating phenomenon first predicted by Maki, Nakagawa, and Sakata , and first applied 
to the solar neutrino problem by Pontecorvo 0. 

The resolution of the solar neutrino problem implies that we now know within certain 
uncertainties the values of Am^2 sin^ 26*12. It is an urgent task to uncover the all possible 
implications of the far more constrained values of these mixing parameters. We address in 
this paper the impact of the KamLAND result on masses and possible mass patterns of 
neutrinos. 

We re-analyze the constraints on neutrino masses and mixing parameters imposed by 
neutrinoless double beta decay experiments in the light of the first results of KamLAND 
experiment under the assumption of Majorana neutrinos. See ref. for the present status 
of the double beta decay experiments. The constraints we will derive 3| are the maximal 
ones that can be derived in a model-independent way (and assuming ignorance of the CP 
violating phases) in a general framework of three-fiavor mixing scheme of neutrinos. It 
utilizes the bounds imposed by neutrinoless double beta decay search and the CHOOZ 
reactor experiment 

A natural question would be "what is the impact of the KamLAND data on the double 
beta decay bound on neutrino mass and mixing parameters?" . The answer is that it allows us 
to derive, for the first time, an upper bound on the neutrino contribution to the cosmological 
Q parameter, the energy density normalized by the critical density in the present universe. 
It is not because the parameter region for the LMA solution becomes narrower (in fact, the 
region for 612 is virtually identical with that before KamLAND), but because the LOW and 
the VAC solutions which allow maximal or almost maximal mixing are eliminated by the 
KamLAND data at 99.95 % CL JL]- 
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II. DERIVING THE JOINT CONSTRAINTS BY DOUBLE BETA DECAY AND 
REACTOR EXPERIMENTS 



In this section we review our formalism presented in [lil . We use the following standard 
parametrization of the MNS matrix 0: 



MNS 



C12C13 S12C13 sise 

-S12C23 - Ci2S23Si3e*'^ C12C23 - Si2S23Si3e*'^ S23C13 



S12S23 - Ci2C23Si3e'" -C12S23 - Si2C23Sl3e"' C23C13 

The mixing matrix for three Major ana neutrinos is 

[/ = f/MNsxdiag(l,e*^e*^), 



(1) 



(2) 



where ^and 7 are extra CP- violating phases which are characteristic to Majorana par- 
ticles 3|. In this parametrization, the observable of double beta decay experiments is 
described as 



(m) 



i=l 

2 2 I 2 2 2i0 I 2 2i(^-5) 



(3) 



where Uet denote the elements in the first row of U and mj = 1, 2, 3) are the neutrino mass 
eigenvalues. In the convention of this paper the normal hierarchy means mi < m2 < ma 
and the inverted hierarchy m^, < nii < 1112. 

In order to utilize (m)™^^ which is an experimental upper bound on {m)pp, we derive a 
theoretical lower bound on (m)^^. With appropriate choice of the phase-factor e^*'-'''"''-' in 
eq. (jni) it leads to the inequality 



{m)Z^ > {m)pp > c 



13 



miCi2 + m2Si2e^*^ 



- "^3Si3. 



(4) 



Strictly speaking, the right-hand-side (RHS) of (jU should be understood as the absolute 
value, but one can show that the alternative case does not lead to useful bound jl^. The 
RHS of Q is minimized by replacing e^*'^ with —1 and sf^ with s^h? respectively, where s^g 



is the largest value of 5^3 which is allowed by the CHOOZ reactor experiment 
speaking, s^h — 0.03. Thus, we obtain 
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Roughly 



(m) 



max \ . 



2 

CH 



miC^2 ~ "^2'5i2 



-CH 



1 1 

— (m2 + mi ) cos 26'i2 - — (m2 



mij 



(5) 



It is convenient to utilize cos 2^12 as a parameter because it arises naturally in the degenerate 
mass case. 

Next, we derive a theoretical upper bound on {m)pp to utilize an experimental lower 
bound {m)™™. Since the RHS of (jHl) is maximized by setting the phase-factors to the unity, 
we obtain 



{m)f^ < (miCi2 + msSis) C13 + m3Si3. 



(6) 



3 



Furthermore, the RHS of (jH)) is maximized by replacing 5^3 with s^h; and zero for the normal 
and the inverted hierarchies, respectively. Then, we obtain 

-{m2 + nil) - ^("^2 - mi) cos26'i2| Cch + "^sSch (7) 



for the normal hierarchy, and 



= y (m2 + mi) - y(m2 - mi) cos26'i2 (8) 

for the inverted hierarchy. 

The obtained bounds (0), and (jHI) can be represented as allowed regions on 

mi - cos 2^12 plane, where mi denotes the lightest neutrino mass for each hierarchy. Notice 
that other mass eigenvalues are expressed by mi together with the solar and the atmospheric 
Am^'s; 



mg = mi + Ami2 + Amgg, 

ml = m\ + Am^2? (9) 



where Am?^ = m^ — m^. 

Alternatively, one can use m/j, the highest neutrino mass (ms for the normal and m2 
for the inverted mass hierarchies) to make the similar plot, as done in 3|. Whichever 
variable is used, they are observables in single beta decay experiments; modulo definitions 
of effective mass parameters [3, 113] it can be related, using ©, the observable in direct mass 
measurement in tritium beta decay experiments. For a review, see e.g., jl^. We emphasize 
that importance of combining informations obtained by the single and the double beta decay 
experiments. 

We do not present a careful treatment of the uncertainty of nuclear matrix elements in 
this letter, but a brief comment will be made on how it affects the mass bound to be obtained 
in the next section. 



III. ANALYZING THE JOINT CONSTRAINTS BY DOUBLE BETA DECAY 
AND REACTOR EXPERIMENTS WITH PRESENTLY AVAILABLE DATA 

Let us first discuss the constraint based on the presently available experimental upper 
bound {m)pp < 0.35 eV (90% CL) [l^. In Fig. □ we display the constraint imposed on 
mixing parameters on m; - cos 2^12 plane, where mi indicates the lowest mass, mi and ma 
for the normal and the inverted mass hierarchies, respectively. The region left to the curves 
is allowed by the inequality (0). The solid and the dashed lines are for the normal and the 
inverted mass hierarchies, respectively. One notices that the two curves completely overlap 
with each other at the present sensitivity of the experiment. The mass squared difference 
parameters are taken as Am^^^^ = 2.5 x 10~^eV^, the best fit value in the atmospheric 
neutrino analysis Amfg = 7.3 x 10~^eV^, the best fit value of the combined 

analysis of the KamLAND and all the solar neutrino data jl8l |. 
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Also shown in Fig. as the shaded strip is the region allowed by joint analysis of the 
KamLAND and all the solar neutrino data [18|- From Fig. ^ one can draw the upper bound 
on mi as mi < 2.13(4.26) eV, where the number in parenthesis here and in eq. ()10|) indicates 
the one for the case where a factor of 2 uncertainty of nuclear matrix elements is taken into 
account. We notice that if either the LOW or the VAC solution were the solution to the 
solar neutrino problem the allowed region of cos 26i2 would extend to zero (maximal mixing) 
and hence there would be no upper bound on mi. 

This bound can be translated into the neutrino contribution to the cosmological Q pa- 
rameter [l^ 

= ^^X' < o.o70(o.i4)/i-2 (10) 

yi.oe V 

which is already very restrictive. Here, h denotes the Hubble parameter in units of 
lOOkm/s Mpc. For possible future use, we provide in Fig. |21 the upper limit of VLy as a 
function of {m)fp^ for three values of tan^6'i2, tan^^^i2 = 0.46, 0.33, and 0.67. These values 
are chosen as the best fit, the upper and the lower limit of 90 % CL allowed region of the 
LMA solution by combined analysis of KamLAND and all the solar neutrino data It 
is notable that even the double beta bound (m)'^^^ < 0.1 eV which may be reachable in 
several years is powerful enough to constrain Qi, parameter to be less than Qi, ~ l/i~^ %. 

We notice that a short (almost) straight line, which is the lower bound on Q^h'^, appears 
at around the lower end in Fig. It is due to the fact that a vertical band appears in 
allowed region in mi - {m)^^ plane at low (m)^^, the feature first observed in |20|. The 
behavior of the curve at around the lower end in Fig. |2b for the case of inverted mass 



hierarchy is slightly different from that obtained in |21| 

We mention here that a similar constraint on neutrino masses may be obtained by cos- 
mological observation 2^ . 



The constraints imposed on neutrino mixing parameters by neutrinoless double beta 
decay have been discussed by many authors. In a previous paper [l^ we tried to classify 
them into the ones in the "early epoch" [i^, the "modern era" i^l, and the "post-modern 
era" l25l which now include some additional references. 



IV. ANALYZING THE JOINT CONSTRAINTS BY DOUBLE BETA DECAY 
AND REACTOR EXPERIMENTS WITH POSSIBLE FUTURE DATA 

In this section, we analyze the constrains based on the double beta decay bound which 
would be obtained in the future. We consider two sample cases. 

Case 1: 0.1eV< {m)pp < 0.3 eV 
Case 2: 0.01 eV < (m)^^ < 0.03 eV. 

Case 1 is a typical mass region proved and to be proved by the ongoing experiments 0, 
0, ll^l, while Case 2 is a typical one for the next-generation experiments. The latter 
include, in alphabetical order, CUORE, EXO, GENIUS, Majorana, MOON, and XMASS 
[28, 29,, ,3^ ,31, 13^, 3^- The mass square difference parameters are fixed as in the previous 
section. 
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A. Case 1 : 0.1 eV < {m)pp < 0.3 eV 

The bounds for this case are presented in Fig. 01 One notices that the lower bound on 
© and (jHl), gives an almost vertical line at nii ^ {'m)^^^ because of the small values 
of Am^2- Therefore, (m)™^"^ is approximately regarded as the lower bound on mi. In Case 
2 with an order of magnitude smaller value of {m)^^ the situation is quite different, as you 
see in the next subsection. 

On the other hand, the upper bound on (m,)^^ gives a nontrivial constraint, as shown in 
Fig. El The curve for (0) approaches to the asymptote cos 26^12 = i^cH large m;, where 
^CH = ■5ch/cch- Hence, if either the LOW or the VAC solution were the case, there would 
be no upper bound on mi. In the LMA solution the upper bound does exist for the current 
allowed region of cos 26*12. 

Notice that the bounds © for the normal and the inverted hierarchies are very similar 
to each other. It means that the degenerate mass approximation holds to a good accuracy. 
In this approximation, the constraint ((H)) becomes 

(m)-- > m [cl^ |cos26^i2| - 4h) • (H) 

where m^ c:^ m. 

We now understand that in the degenerate mass regime mi/ (m)™^^ is a good parameter. 
In Fig. m the bounds © are presented on the mi/{m)'^p^ - cos 26*12 plane for {m)'^^^ = 
0.1 eV. Two curves for the normal and the inverted mass hierarchies overlap well, indicating 
goodness of the degenerate mass approximation. Once we know that the scaling works, the 
result presented in Fig. jH applies to any values of mi and {'m)'^p^ as far as the degenerate 
mass approximation holds. 



B. Case 2 : 0.01 eV < {m)pp < 0.03 eV 

The bounds for this case are presented in Fig. jSl It is clear that the degenerate mass 
approximation is no longer good because the relevant energy scale is smaller than the at- 
mospheric mass scale y^Am^m — 0.05 eV.^ The bounds for the normal and the inverted 
hierarchy significantly differ with each other. 

An important point in Fig. El is that with such improved sensitivities experiments may 
start distinguishing the two different mass hierarchies. Most notably, the difference between 
the normal and the inverted mass hierarchies is quite visible at m/ < 0.05 eV, and the lower 
limit of m; disappears for the inverted hierarchy. It occurs because the constraint (|HI) is 
satisfied even for mi = 0. A similar things occurs also in the normal hierarchy if (m)™™ is 
sufficiently small. One can easily estimate, by taking mi = and cos 26*12 = in and 
the critical value of {m)™™ under which the lower bound on mi disappears. Approximately, 
they read 

^ c^jj^A^ + s^HV^A^atm ^ 0.0056 eV ~ (12) 



Although we put aside the uncertainties of Am^'s for simphcity, those affect particularly the values 
related to the lower bound on m; in such a small mass scale. Therefore, the values in (|12|l . (|13|) . and the 
left-hand-side of \ib\ are just for reference. 
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for the normal hierarchy, and 

4h (^^AmL. - Amf2 + ^j ^ml^ ~ 0.048 eV ~ ^^mL, (13) 

for the inverted hierarchy. In Case 2, (m)™^"^ = 0.01 eV is smaller than 0.048 eV, and that is 
why there is no lower bound on m; for the inverted hierarchy in Fig. El 

Finally, let us extract the bound on rrii from the result of js^] by assuming that it will be 
confirmed by NEMOS. Since {m)^p^ = 0.84 eV is large enough, we can use Fig.HJ Therefore, 
the upper bound on m/ is extracted as 

mi < 3[6] X {m)f^'' = 2.5 [5.0] eV (14) 

for the LMA best fit parameters and for the most conservative case in square parenthesis, 
respectively. On the other hand, {m)™^^ = 0.05 eV determines the lower bound on m;. By 
combining those results, we obtain 

0.05 eV(0.005 eV) < mi < 2.5[5.0] eV (15) 

for the normal (inverted) hierarchy with, in right- hand- side, the LMA best fit parameters 
and for the most conservative case in square parenthesis. Therefore, we see that all the 
active neutrinos is massive if the result of 34] is confirmed. This conclusion, however, may 
not survive if the uncertainties of Am^'s and nuclear matrix elements are taken into account. 

V. SUMMARY 

In the light of the KamLAND data we reanalyzed the constraint on neutrino masses 
and mixing imposed by neutrinoless double beta decay search and the CHOOZ reactor 
experiment. We have shown that by now the absolute mass scale for neutrinos is bounded 
from above by the "conspiracy" of KamLAND and double beta decay experiments. This 
conclusion and our whole analysis in this paper are valid under the assumption that neutrinos 
are Majorana particles. 

In closer detail, allowed regions on a plane spanned by the lowest neutrino mass versus 
the solar mixing angle 612 were obtained by using (m)™|'^ ((m)™^"^), the experimental upper 
(lower) bound obtained and to be obtained in the ongoing and the future double beta 
decay experiments. For given 612, these become constraints on neutrino mass, the lightest 
mass mi (or the heaviest one though we did not give them explicitly). Roughly speaking, 
(m)™™ < TTii < 3[6] X (m)™^^ for the LMA best fit parameters and for the most conservative 
case in square parenthesis, respectively. 

It was made clear that the condition | cos 26^12 1 > — 0-03 was necessary for the upper 
bound on mi to exist. On the other hand, the condition (m)™™ > 0.0056 eV(0. 048 eV) needs 
to be satisfied for the normal (inverted) mass hierarchy, for the lower bound on mi to exist. 
For example, 0.05eV < {m)i3p < 0.84eV gives 0.05 eV(0.005 eV) < < 2.5[5.0]eV for the 
normal (inverted) hierarchy. 

Note added: 

While this paper was in the revision process, we became aware of the report by the 
Wilkinson Microwave Anisotropy Probe (WMAP) j3^ in which a severe constraint on omega 
parameter Q^h^ < 0.0076 (95 % CL) is placed. It means that double beta decay searches 
must have the sensitivity in the region {m)^^ < 0.1 eV to be competitive, as one can see from 
Fig. 2. 
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°0 0.5 1 1.5 2 2.5 3 

m/(eV) 

FIG. 1: The double beta decay constraint with the CHOOZ bound (jSJ is displayed as an allowed 
region in rrii - cos 2^12 space. The region left to the curve is allowed. The solid and the dashed lines 
are for the normal and the inverted mass hierarchies, but they completely overlap. The allowed 
LMA parameter region by the combined solar-KamLAND analysis is superimposed with shadow, 
and cos 2^12 = 0.37 is the best fit value. 
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FIG. 2: Plotted is the upper bound on $1,^, the neutrino contribution to cosmological J7 parameter. 
Three hnes correspond to, from above to below, tan^ 9i2 = 0.67, 0.46, and 0.33, respectively. Fig. [2^ 
andl^Ja are for the normal and the inverted mass hierarchies, respectively; These correspond to the 
edges of the LMA region and the best fit value. The bold horizontal line at Ujyh? = 0.072 is the 
upper bound determined by the recent result of single beta decay experiments. The thin horizontal 
line is the lower bound obtained by mi = 0. 
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FIG. 3: The double beta decay constraints with the CHOOZ bound ©, and © displayed for 
0.1 eV < {m)j3f3 < 0.3 eV. The solid and the dashed lines are for the normal and the inverted mass 
hierarchies, respectively. The allowed LMA parameter region by the combined solar-KamLAND 
analysis is superimposed with shadow, and cos 2^12 = 0.37 is the best fit value. 




mi/< m 



FIG. 4: The double beta decay constraint with the CHOOZ bound © is displayed for (m)™^^ = 
0.1 eV to indicate the scaling behavior of the bound with variable mi/{m)^p^ as used for abscissa. 
The region left to the curve is allowed. The solid and the dashed lines are for the normal and 
the inverted mass hierarchies, respectively. The allowed LMA parameter region by the combined 
solar-KamLAND analysis is superimposed with shadow, and cos20i2 = 0.37 is the best fit value. 
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FIG. 5: The double beta decay constraints with the CHOOZ bound ©, and © displayed for 
0.01 eV < {m)pp < 0.03 eV. The solid and the dashed lines are for the normal and the inverted mass 
hierarchies, respectively. The allowed LMA parameter region by the combined solar-KamLAND 
analysis is superimposed with shadow, and cos20i2 = 0.37 is the best fit value. 
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